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MOMENTS  FOR  A MODIFIED  MAKEHAM  LAW  QF  MORTALITY 


R.  Clifton  Railey,  Ph.D. 

1.  INTRODUCTION 

The  Makeham  model  for  the  force  of  mortality  has  intuitive  appeal  and 
an  excellent  reputation  as  an  appropriate  model  for  use  in  actuarial  problems 
(Jordan,  [1967]).  Recently,  I have  been  interested  in  survival  problems 
associated  with  medical  procedures  for  which  the  force  of  mortality  following 
the  procedure  is  initially  high  and  eventually  falls  to  a lower  rate.  Since 
follow-up  is  generally  limited,  long-term  extrapolation,  which  may  involve  an 
increasing  force  of  mortality  at  some  longer  period,  is  not  modeled  beyond 
the  constant  term. 

I offer 

h(t)  = 6 + aexp(-yt)  (1) 

as  the  force  of  mortality  for  what  I call  the  modified  Makeham  model.  Note 
that  the  parameter  6 is  the  long-term  risk,  c*  is  an  initial  excess  risk,  and 
y is  a rate  constant  for  the  disappearance  of  the  initial  excess  risk.  For 
this  force  of  mortality  the  probability  of  survival  is 

l-F(t)  = exp[-6t-(oi/Y)  (l-exp(-yt))  ] . (2) 

In  some  trial  analyses  maximum  likelihood  estimation  has  been  used  to 
obtain  estimates  of  the  parameters  for  this  model.  The  likelihood  function  was 
set  up  to  reflect  that  not  all  patients  fail  in  the  course  of  the  study.  The 
iterative  maximum  likelihood  procedure  we  used  employs  a routine  for  numerical 
derivatives.  This  has  facilitated  the  exploration  of  alternative  models.  The 
following  is  a development  of  a computationally  useful  expression  for  computing 
the  moments  for  this  model. 
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2.  MOMENTS 


It  is  shown  in  Appendix  I that  the  v-th  moment  about  the  origin  is 
given  by  , 

H'  = vj-  fcV"  [i-F(t) ]dt,  v = 1,2,...  (3) 

v 0 

for  a non-negative  random  variable  with  cumulative  distribution 
function  F(t)  provided 

lim  tV(l-F(t) ) = 0.  (4) 

For  the  modified  Makeham  model  the  required  limit  is  shown  in  Appendix  II. 

Now  these  results  can  be  applied  to  the  problem  of  finding  the  moments 
about  the  origin. 

That  is,  m ^ 

p'  = v/  tV  exp[-6t-(a/y) (l-exp(-yt)) ]dt.  (5) 

v 0 


Next  use  the  exponential  expansion  for  exp[  (a/y)exp(-yt) ] and  do  the 
integration  term  by  term.  This  gives 

u"  = vexp(-a/y)  £ J°°  tV_1exp[ -(6+iyt]dt 

i=0  0 

= vexp(-a/y)£  — ^ — 

i=0  (6+iy) v 

CO 

_ r(v  +1)  exp(-a/Y)  l (a/y) 1/[  il (1+iy/ 6 ) V ] , (6) 

6V  i=0 


where  the  gamma  function  T(v+1)=  v!  for  v = 1,2,3,.... 

The  raw  moments  given  in  equation  (6)  are  used  to  compute  the  central 
moments  in  the  usual  way. 

That  is. 


u 


V 


p 


v-i 
1 ‘ 


(7) 


2 


Prom  the  central  moments  the  skewness 


- V> 3 / y 2 


and  the  kurtosis 


82  = ^4 /y J 


are  computed. 


3.  EXAMPLE 


In  a trial  data  set  the  maximum  likelihood  estimates  for  a,  y,  and  6 wer 

<\  * a 

a » 1.388,  y = 3.506,  and  6 = 0.120.  The  expressions  given  above  were  used 
to  compute  estimates  of  the  first  four  moments,  the  central  moments,  the 
skewness,  and  the  kurtosis.  The  results  are  for  the  moments 


yf  = 5.691, 


y 2 = 93.526, 


y 3 = 2,337.114, 


U4  = 77,902.685; 


for  the  central  moments 


y 2 = 61.143, 


y 3 = 1109.016, 


P4  = 39,730.361; 


for  the  skewness 


- 2.320; 


and  for  the  kurtosis 


B2  - 10.627. 
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The  calculations  for  the  moments  were  done  by  Mr.  James  P.  Summe  on  a 


I I 


FORTRAN  program  which  he  developed. 
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APPENDIX  I 


Let  p'  be  the  v-th  raw  moment  about  the  origin.  Then  some  elementary 
operations  can  be  used  to  verify  that 

u'  = vf  tV  1 [ 1-  F(t)  ]dt , v = 1,2,... 

v 0 

for  a random  variable  t with  cumulant  distribution  fuction  F(t)  . This 

may  be  verified  directly  through  integration  by  parts  if  we  put  u-  l-F(t) 

and  dv  = tV  Fdt  and  if  F(t)  is  such  that 

lim  (l-F(t))tV  = 0. 
t-*® 

The  result  also  follows  from  a development  based  on  the  moment-generating 
function 

M (s)  = f°  eSt  h(t)exp[-ft  h(x)dx]dt. 

0 0 


ve  set  u=  eSt  and  dv  = h(t)exp-Jt  h(x)dx. 


re  s is  such  that 

5t  ft 


lim  e exp-/  h(x)dx  = 0. 


t-*“ 


That  is,  we  require 


s </*"  h(x)dx 
0 


for  some  t.  Then 


M(s)  = 1 + s/°°eStexp[ -/fc  h(x)dx]dt. 
0 0 


From  this  we  find 


dV  M(s)  - v f tV-1eSt(i-F(t»  + sl0  t e (l'F(t))dt 


or  for  s-0. 


d s 


u'  = tv  1[l-F(t)]dt  V = 1,2,.., 
v 0 


In  the  same  way  the  central  moment-generating  function  can  be  used 

to  show  that  the  v-th  central  moment  is 

, v— 1 . 


g * vj  (t-m)V-i(l-F(t))dt. 


( 

! 
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APPENDIX  II 


The  theorem  for  the  moments  can  be  applied  to  our  problem  if  we  can 
show  that 

L = lim  (N(t)/D(t))  = 0, 

t-KO 


where 

N(t)  = t\ 


and  D(t)  = Q(t)exp(6t), 

when 

Q(t)  = exp[  (a/y) (l-exp(-yt)) ]. 

Note  that 

N(v)(t)  = v! 

for  v = 1,2,....  Also  observe  that  Q(<=°)  = exp(a/y)  and  that  by  induction 
we  can  show  that  (°°)  * 0 for  i = 1,2,.... 

To  show  that  Q^(°°/  * 0,  we  first  note  that 
Q'(t)  - aQ(t)exp(-yt) . 

Rrom  this  we  obtain 

q(1+1> (t)  - aexp(-yt)  S1(t), 


where 


t(t)  - l Q(j)(t)(-y)i-j. 

1 j-0  ^ 


If  we  assume  Q^(“>)  * 0 for  i>l,  then 


Sj^®)  - a(-y)  exp(ct/y) , 


which  Implies 


Q<1+1)(“)  - 0, 


since  is  finite  and  exp(-fit)-K)  as  t-*». 
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Now  a similar  argument  applied  to  D(t)  gives 


where 


frcai  which  we  have 

lim  R^Ct)  = S^expCa/y). 

Consequently, 

lira  D(i) (t)  = ®,  i - 0,1,2 

t-KK> 

Repeated  application  of  L'Hospital's  rule  gives 


L 


lim  N(v) (t)  - 0 
D(v)(t) 


when  v is  an  integer. 


